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^pr ■ Abstract 

_~ , A model of the mean-field interacting boson gas trapped by a weak harmonic 

j^* potential is considered by the boson random point fields methods. We prove that 

^^ ■ in the Weak Harmonic Trap (WHT) limit there are two phases distinguished by 

the boson condensation and by a different behaviour of the local particle density. 
For chemical potentials less than a certain critical value, the resulting Random 
Point Field (RPF) coincides with the usual boson RPF, which corresponds to a 
S^ . non- interacting [ideal) boson gas. For the chemical potentials greater than the 

critical value, the boson RPF describes a divergent (local) density, which is due 
to localization of the macroscopic number of condensed particles. Notice that it is 
this kind of transition that observed in experiments producing the Bose-Einstein 
Condensation in traps. 
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1 Introduction and Main Results 

1.1 Weak Harmonic Traps 

We consider the quantum statistical mechanical models of boson gas equipped with a 
K-parameterized family of one-particle Hamiltonians of harmonic oscillators: 

iA/ a 2 ^ i\ 



K 2 ^ V dx* K 2 

which are self-adjoint operators in the Hilbert space Sj := L 2 (M. d ). 

In this setup a "thermodynamic limit" corresponds to k —> oo (i.e. the "opening" 
of the trap [DGPS]), which we call the Weak Harmonic Trap (WHT) limit. Notice that 
the set C^°(M. d ) is a form-core of the operator (1.1) and that this set is also a form-core 
for the operator (— A)/2. Here A is the standard Laplace operator in M d . Then (see 
e.g. [Ka]) one obtains the strong resolvent convergence: 

lim h K = (-A)/2 . (1.2) 



In spite of convergence (1.2), there is a drastic difference between the properties of the 
infinite Ideal Boson Gas (IBG) systems prepared via the WHT limit and via standard 
thermodynamic limit (TDL) liin^oo A^ = M. d with the one-particle operators ti '■ = 
(— A/2)x, with a "non-sticky" (e.g. Dirichlet) boundary conditions [VVZ]. Here A^=i C 
M. d is a bounded open region of unit volume |A^ = i| = 1 containing the origin x = 
whose boundary 8Al=i is piecewise continuously differentiable and 

A L := {x E R^L^x E A L=1 } , L > 0. (1.3) 

In fact, it is known that the Bose-Einstein Condensation (BEC) occurs for dimensions 
d > 1 in the IBG via WHT limit k — > oo, while for dimensions d > 2 in the IBG via TDL, 
see (1.13) and [DGPS], [PeSm]. Similarly, it is well-known that thermodynamic proper- 
ties of the boson gases are very sensible to different ways of taking the thermodynamic 
limit [LePu] or to the different choices of the boundary conditions [VVZ]. 

The purpose of the paper is to examine the position distribution of the mean field 
boson model in WHT limit and to compair its behaviour to those of the mean field 
boson models in TDL or ideal boson models in WHT limit. 

Our method is based on the theory of Random Point Fields (RPFs) (see e.g. [DV]). 
The usual boson and the fermion RPFs [Ly, M75, M77] have been formulated in a 
unified way in terms of the Fredholm determinant together with other related RPFs, 
which are indexed by fractional numbers, in [ShTa]. They have been re-derived as 
theories which describe position distributions of the constituent particles of quantum 
gases in the thermodynamic limit for canonical ensembles in [TIa]. It was shown that 
the Random Point Field (RPF) corresponding to fractional numbers [ShTa], describes 
the gases which consist of particles obeying the parastatistics [Tic] . The RPF describing 
a homogeneous Bose-Einstein condensation have been studied for the first time in [Tib], 
where the RPF is given by the convolution of usual boson RPF and another RPF. The 
latter one seems to describe position distribution of the condensed part of the constituent 
bosons. This RPF has been re-formulated using the Cox process. [EK] 

These theories of RPFs yield a precise information about the position distribution 
of the constituent quantum particles, although they are not suitable to characterize the 
quantum systems completely (however see [FF, F]). 

1.2 Ideal and Mean-Field Boson Gases in the WHT 

Recall that in the grand- canonical Gibbs ensemble the partition function of the IBG 
trapped by harmonic potential (1.1) is given by 

oo 

E 0tK (/3,fi) := J2^ n Tr^ mm [^G K (/3)] . (1.4) 



n=0 



Here f)™ ymm : = (<& n L 2 (R. )) synirn is the n-fold symmetric Hilbert space tensor product of 
JO := L 2 (R d ), G K (f3) = e~ /3h,K the one particle Gibbs semigroup. The zeroth term in (1.4) 
equals to 1 by definition. We consider the case of positive inverse temperature j3 > 
and of negative chemical potential \x < 0. 



The spectrum of the operator (1.1) is discrete and has the form: 

Spec(/i K ) = { e K (s) := \s\i/k \ s = (si, • • • , s d ) G Z+ } (1.5) 

where |s|i := Y^i=i s j anc ^ ^+ ^ s the se ^ °f an non-negative integers. The normalized 
eigenfunctions of the states for s G Z+ are given by 

1 d 1 

0«,«(a:) = -rf/i </>,(x/</k) = Yl -T/1<Ps 3 {xj/\/k) , (1.6) 

i=i 
where for each component j, <p s . is related to the Hermite polynomials H s .(z) by 

<P S] (x,) = (2 s * Sj \ ^r 1/2 H S] ( Xj )e-^/\ j = 1, . . . , d . (1.7) 

The ground state is denoted in this paper by 

n «G0 = r^wi e " N2/2K = <f>*=oAx) , (i.8) 

where x = (xi, ■ • • , a;^) G M d , |x| 2 := X/j=i x i- 

Integral kernel of G K (j3) = e~ l3tlK has the explicit form (the Mehler's formula for 
oscillator processes): 

exp{-(2/ t )" 1 tanh(/3/2/ t )(|x| 2 + \y\ 2 ) - \x - y\ 2 /(2 K smh(f3/K))} 
^KP^V)- {tt«(1 - e-W«)}d/2 • U- y J 

Here the operator G K ((3) belongs to the trace-class £i(L 2 (IR d )), with the trace-norm 
equals to TrG K ({3) = 1/(1 — e~ /3 ^ K ) d = 0(n d ) for large k. The largest eigenvalue of 
G K ((3) coincides with the operator norm ||G K (/?)|| = 1. We write all the eigenvalues of 
operator G K (j3) in decreasing order: 

9 { o K) = 1 > 9 [ K) = e-^ > g { 2 K) >■■■. 
The expectation value of total number of particles is given by 

N , M .1^^^_±_. (1 , 0) 

Since the value (1.10) diverges in the WHT limit k — > oo as n d , one introduces the 
scaled quantity [DGPS], [PeSm]: 






which is a Darboux-Riemann sum for the integral 

p(/5,/i)= limp K (/5,/i)= /" -^ 





- 1 
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Since n d may be interpreted as the effective "volume" {cf. Remark 1.1), p(/3,p) and 
p K ((3, p) are regarded as the expectation value of effective space-averaged density of the 
system, which has non-homogeneous space distribution. One defines its critical value as 
usual: 

p c (P):=su V p(P,p)=((d)/p d . (1.13) 

Notice that (1.13) is bounded for d > 1. Therefore, if p > p c (P), the IBG in the 
WHT limit k — ► oo manifests a BEC in the scaled oscillator ground state (1.8) with the 
expected space-averaged condensate density: 

PM ■= P - pM = Hm ^ d e . r J M _ 1 • (1.U) 

Moreover the expected local density po(/3)(x) can be defined such that 

ft<fl(»> = is ? «-£<£!-! (1 - 15) 

holds. Here ~fi K ((3,p) is the unique root of the equation, p = p K ((3,p) cf. (1.11). Their 
limits satisfy 

KmJI K {P,p<p c (P))<0 and lim /Z K (/3,p ^ p c (/?)) = 0, (1.16) 

especially 

R (ftp> feW ) = - /?(/) _^ )K + o( K -). (1.17) 

The integrated density of states N K (E) for the operator h K is given by 



rv 



Then from its Laplace transform 

i 

e~ tE dN K (E) = [k(1 - exp(-t/K))]" 



'o 
we obtain the k — » oo limit 



d/V(£) = ^ dE . (l.is) 



In terms of these density of states, (1.11) and (1.12) are written as 



oo 



dM K {E) f°° dAT(E) 



It is instructive to compare these results with properties of the IBG "prepared" via 
standard thermodynamic limit L — > oo (1.3) for e.g. Dirichlet boundary conditions. It 
is well-known [ZB] that in this case the expected boson density is 

PaMu) = TO E e ^fa-)-M) - i • (L19) 



Here {£l(J)}j^o is the spectrum of the one-particle operator t L and /i < min^ £l(J) ~^ 
as L — > oo. It also has the expression 

f>AM»)=£ X-f-i (1 - 20) 

in terms of the integrated density of states {Nl{E)}l >q . The thermodynamic limit 
J\f(E) = limi^oo Af L (E) is independent of "non-sticky" boundary conditions [RSIV] 
and given by 

E (d-2)/2 

Ml - E) = Wr W 2) dE - (L21) 

Then we get the limit of the expected density 

mri = gnPAMn) = [ J^.ifl i ■ (i-22) 

The critical particle density for the IBG is 

&(/?):= sup p(/3 )A i)= r^\ = C(d/2)/(2n^. (1.23) 

m<o Jo & - 1 

Note the difference between p c (j3) and p c ((3). In particular (1.23) is bounded only when 
d > 2. Thus thermodynamic properties of the IBG in the standard TDL A^ — > M. d and 
the WHT limit k — > oo are different in spite of delusive impression that they have to 
produce identical systems. 

Now we consider the mean-field interacting bosons trapped in the harmonic potential 
(1.1). Its grand- canonical partition function is given by 

oo 

S^GM := J2e^ n ~ Xn2/2Kd) Tr^ ym J® n GM] . (1.24) 

n=0 

We consider the case of (3 > 0, A > and arbitrary /j, G M. Hereafter, we suppress the 
symbol A from the left-hand side of (1.24), since we fix A > in the rest of the paper. 

Remark 1.1 The scaling with the "volume" n d (imposed by the WHT) is a conven- 
tional way to consider the Bose-Einstein condensation in traps, see [DGPS], [LSSY], 
[PeSm]. Our definition of the MF interaction in WHT applies a space-average over the 
"volume" n d , which plays the same role as \Al\ in the standard mean field model where 
the interaction has the form An 2 /2|Ai|, see e.g. [ZB]. 

Notice that A > corresponds to repulsive mean-field (MF) particle interaction, 
whereas A = is the case of the IBG (1-4). 

In the present paper we consider in (1-24) only the "square" mean-field repulsive 
interaction U<i := An 2 /2n d . Although application of the Large Deviation technique makes 
it possible to consider also the case of the general MF interaction U<& := K d <&(n/ n d ), 
where $ : R — > K is a piece-wise differentiable continuous function bounded from below, 
see [TZ]. 
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To study the non-homogeneous condensation and the space distribution of the con- 
stituent bosons in the system (1.24) we use the RPF z/ K/ g, i.e., the probability measure 
on the space of locally finite point measures with generating functional: 

iW/) : = ^A e ~ m ] ( L25 ) 

= zMri ^ ^ n ~ PXn2 ' 2Kd ^ *— i® n (G^)e-f)} , 

where / G C*o(lR d ), / ^ 0. Here E Kii g i/i [ • ] stands for expectation with respect to 
v K ,/3 t n(d£), and £ denotes the integral variable which represents locally finite point mea- 
sure, see [DV, TIa, Tib, Tic]. The measure v K $^ describes a finite RPF whose Janossy 
measure can be given explicitly, see Remark 2.1. 

In the present paper we study the properties of the MF interacting boson RPF v K ,p,n 
in the WHT limit k — ► oo by analyzing the generating functional (1.25). To this end 
we first define the MF critical chemical potential 

, cm \f dp , (°° dM(E) AC(rf) ,,,-> 

This critical parameter is similar to the critical chemical potential Jl c (l3) := Xp c (/3) for 
the standard homogeneous MF boson gas via TDL, see e.g. [ZB] and (1.23). 

1.3 Main Results 

Now we can formulate our main theorem: 

Theorem 1.2 (i) Let fi < fi\ tC (l3) (normal phase). Then the RPF v k ,p,ij, defined by 
(1.25), converges weakly in the WHT limit k — > oo to the boson RPF vp, r * corresponding 
to the generating functional: 

Ep,rSf)- = ®/3,r*[e- {m ] (1.27) 

= Det [1 + y/l - eS r„ G({3) (1 - r* G^))" 1 y/l - e~f ] _1 . 

Here E^, [ ■ ] denotes expectation with respect to the measure vp jTt , Det stands for the 
Fredholm determinant, G((3) := e^/ 2 is the heat semigroup on S) and r* = r*(/3,/i, A) G 
(0, 1) is a unique solution of the equation : 

ft.-logr + AgjT ^'ffi^ . (1.28) 

(ii) For /j, > fi\ tC ((3) (condensed phase) , the generating functional (1.25) has the following 
asymptotics: 

^^ogE Ki0 M) = ~ ~ ~^ (V^^,(l + K f rVT^ I f), (1.29) 

where K f := (g((3) x ' 2 (1 - G{P))~ l l 2 y/l- e~fX (gW /2 (1 - G((3)Y l l 2 ^l - e~A is 
a positive trace- class operator on fj = L 2 (WL d ) for d > 2. 
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Remark 1.3 For dimensions d > 1 the integral (1.26) is finite. This gives an idea 
that the BEC is possible for the case fi > /^\ )C (/3) and dimensions d > 1. However, in 
the present paper we assume d > 2 to be able to prove our main theorem. In fact for 
d = 2, the operator Kf is not a trace- class operator and the Fredholm determinant is 
not well-defined, see [Tib]. 

Remark 1.4 Because of a technical difficulty, actually we do not have results for the 
critical case \i = fi\ <c (j3). 

Remark 1.5 The heat semigroup G(j3) is appeared here since its kernel 

G(f3;x,y) = (2<K{3)- d/2 e-\ x - y \ 2/2 P 

is the point-wise limit of the Mehler kernel (1.9), as k — > oo. Therefore the generating 
functional of the resulting RPF (1.27) in (i) has exactly the same form as that for the 
standard homogeneous IBG in the non-condensed phase, see eq.(2.13) in [TIa], see also 
[ShTa]. That is, the position distribution of the model in the non-condensed phase is 
coincides with that of the standard IBG. However, it is the integrated density of states 
M(E) that is appeared in (1.28), instead of M(E). It implies that the dependence of r* 
as a function of parameters, especially (3, is different from standard IBG or the standard 
MF-model. 

Remark 1.6 The RPF describing BEC for homogeneous IBG is given by the convolu- 
tion of two RPFs [Tib] . One convolution component is the usual boson RPF, while the 
other component seems to describe the position distribution of the condensed part of the 
constituent bosons. The behaviour of the generating functional in case (ii) shows that 
the latter component overwhelms the former in the present model. It is to be noted here 
that the latter can not be explained by the particles in the ground state alone, it contains 
the effect of the interference between "the condensed part" and "the normal part" [EK], 
although the intensity of the RPF is proportional to the square of the ground state wave 
function, as (1.15). 

The sharp contrast between two regimes (i) and (ii) in Theorem 1.2 may be seen by 
the expectation values. 

Corollary 1.7 For the case (i) \x < p\ tC (/3) (normal phase ) 

E ftr ,[(/,0]=Tr[/r,G(/3)(l-r,G(/3))- 1 ]= Pr , / dx f(x) 

Jm d 

holds, where p rt is given by 

oo 

Pr. = r.G(J3)(l - nG{(3))-\x,x) = ^/^r^ 



d=i 



For the case (ii) /i > p\ tC (fl) (condensed phase), 



'^ ''"f'^ ^lr 1 / *"(«> d- 3 ") 



holds, where f G C ( 



The weak limits of the RPFs concerns the limit of the local position distribution 
of particles. In this sense, the results of Theorem 1.2 and Corollary 1.7 in regime (i) 
may be interpreted as follows: in the WHT limit the position distribution of the MF 
interacting bosons in neighbourhoods of the origin of coordinates (i.e. the bottom of 
the WHT potential) is close to that of a free IBG corresponding to the unconventional 
parameter (1.28). The information about the particle position distribution in domains 
distant from the bottom of the WHT are missing in the limit vp, r * ■ in order to take this 
"tail" particles into account we use the standard definition of the grand-canonical total 
number of particles for our model : 

(tot), a x 1 dlnE K {P, ft) 
PV(/M == ^ ^ (1-31) 



1 oo 

— — — V n e ^n-XnV2^) Tr ^n G 



Since n d is interpreted as the effective volume of the model, p K ° x (j3,fj,) represents the 
effective total space- averaged density of the non-homogeneous system (1.24). 

Theorem 1.8 The WHT limit of (1.31) 

pf 4) (/3,/i) = lim p { *f(P,p) = lim «T d Tr [ n G K (l - nG^' 1 } (1.32) 

exists and satisfies : 
(i) (// < Hx,M) 



p[ tot \(3,p) = [ j^P - and (3p = logr, + \f3pf* } (/5, p) ; (1.33) 

(ii) (p > fJt\M) 



p[ tot) (P,p)=p/\. (1.34) 

It also holds that 

p? ot) ((3):= lim . p[ tot) GM = C(d)//3 d . (1.35) 

Remark 1.9 The readers should not to confuse two "densities" : p Tt in Corollary 1.7 
and p x ° ((3,p) defined above. The p Ttt can be interpreted as the limit of the "local" 
density around the origin of coordinate of non-homogeneous RPF v R $^, on the other 
hand the p\° (fi,p) retains the information about the expectation of the total number of 
particles with respect to v K fi^ through the WHT limit. 

Note that p { c° t] ((3) coincides with (1.13) of IBG in the WHT. 

Remark 1.10 Qualitatively different behaviour of the space distributions of bosons de- 
scribed in Theorem 1.2 can be understand heuristically with the help of Theorem 1.8 in 
the following way. Consider the WHT limit k — > oo; 
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In case (i) the bosons are distributed almost uniformly in the region of radius k 
according to the kernel (1.9). 

On the other hand, in case (ii) [condensed phase) the condensed part of particles 
K d (p x ° (P,fj) — Px° c (P)) = K d {p — px :C (f3))/X is localized in the region of radius 0(k 1 ^ 2 ) 
according to the profile of the square of the ground state wave function Q K . Whereas the 
particles outside of the condensate essentially spread out over the region of radius k. 

The paper is organized as follows. Preliminary estimates and results concerning the 
WHT limit for the mean-field interacting boson gas (A > 0) are collected in Section 2. 
Section 3 and 4 are dedicated to the proof of Theorem 1.2 and 1.8, respectively. We 
reserved Section 5 for summary and conjectures. 

2 Preliminary Arguments and Estimates 

In this section, we write the expectation (1.25) as the ratio S K (/3, p)/E K (/3, p). The 
representations of S K (/3, p) and H K (/3, p) are given in the form of integration of Fredholm 
determinants. We also give the miscellaneous estimates needed for the evaluation of 
these integrals. 

2.1 H K (/?,/i) and E K (j3,n) 

In terms of the projection operator on fj n = <& n L 2 (R d ) onto its subspace jj™ smra , the 
grand-canonical partition function can be written as 

Z K (P,p) = ^e^-^l^- J2 Tr® nL i m [(® n GM)U(a)], 

n=0 <76<S n 

where the second sum is taken over the symmetric group S n and 

U(<x)<pi ®---®(p n = <Ar-i(i) <S> ■ ■ ■ <S> <p ff -i(n) for o- E S n , ip h ■ ■ ■ , ip n G L 2 (R d ) . 

Hence we have 

E K ([3,p) = y^ j Y] / {\\G K (f3)(x j ,x a -i {j) ))dxi---dx n 

n=o n! ^Jm n f=i 

= Yl ~\ / ^ { G M( X ^ X j)}l^,^n dXi ■■■ dx n, 



n=0 



here Per stands for the permanent of the matrix {G k (/3)(xj,Xj)} . 

Remark 2.1 The point field Vk^^ of (1.25) can also be defined in terms of Janossy 
measures or exclusion probability [DV] . This means that v K p ^ is a finite point field, 
which assigns the probability 



e /3nn-f3\n 2 /2K d 



Pr{dX n } := — ,_, — - Per {G K (/3)(xi, x j )} 1 ^ J<n dxi ■■■dx n 
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to the event {dX n }: there are exactly n points, one in each infinitesimal region 

As in [TIa, Tib], we use the generalized Vere- Jones' formula [ShTa, V] in the form 

— / PerfJ(xi,Xi)\i<ij< n dx\ ■ ■ -dx n — ® — — — : -, 

n\ J XK jJf^hj^ri n J Sr{0) 27r iz n+1 Det(l - zJ)' 

where r > satisfies \\rJ\\ < 1. S r (Q denotes the integration contour defined by the 
map 9 *— > £ + r exp(i9), where 9 ranges from — it to it, r > and ( G C. Then we obtain 

"^ h 1(0) 2^^Det(l - zGMY 

where r G (0, ||G f K (/3)|| _1 ) = (0, 1). Note that the zeroth term is 1 in this expression. 
Let us substitute 



-l3\n 2 /2K d _ 


/ px 

V 1lTK d , 


// xexp (-i^ ( 


If s > satisfies 




e f3/i-f3\s/K d < r 



■j((x + is) 2 - 2in(x + is))) . (2.1) 



we can take the summation over n together with the complex integration and a scaling 
of x to get 



2vTK d 7 E J 2iri [z - e ^+^^- s )/ Kd )Det[l - zG K ((3)} 



K 



d (■ -isx-K d x 2 /2(3\ 



e^l^ d / dx r \ . —— d — - . (2.2) 



2np\ J R Det[l-e^+ ix -^ s / Kd G K (/?)]' 

Note that after ^-integration, r disappears and (2.2) is valid for any s satisfying exp(/3/i — 
(5\s/K d ) G (0, ||G K (/3)|| _1 ). We will estimate the integral in the spirit of saddle point 
method. Here, we extract the main part from the integral. Let s = s K ,r = r K be the 
unique solution of the system: 

r = exp (/?// - p\s/n d ) 

s = Tr[rG K ([3)(l-rGM)- 1 }. 

Obviously, the condition r K G (0, ||G K (/3)|| _1 ) is fulfilled. Hence, we can substitute in 
(2.2) s by s K . Using the product property of the Fredholm determinant, we get for 
denominator of (1.25) the representation: 



K d e (3\sl/2 K d r e -is K x-K d x 2 /2l3\ 

£ K (/3,/i) = «/ - QX - / dx 



2ir(3\ Det[l - r K G K ((3)] J R Det[l - (e« - l)r K G K (/3)(l - r K G K {P))-i\ 

(2.4) 



11 



For the numerator of (1.25), we introduce bounded symmetric operators 

GM(f) ■= GM 1/2 e~ f GM 1/2 , (2-5) 

indexed by function / G Co(K d ), / ^ 0, which we skip below for simplicity Then for 
generating functional (1.25) one gets the form: E K ^^(f) = E K ([3, /i)(f)/E K (/3, fi), where 

oo 

s«0M(/) = £«6M = Y. e ^ n ~ pXn2/2KdTT ^J® n GM] 

n=0 



K d p f3Xs 2 K /2K d f p -is K x-K d x 2 /2f3\ 

- / dx - - - (2 6) 

2tt/5A Det[l - f K G K {/3)] A Det [! " (^ ~ 1)?AM(1 - W))" 1 ] 
Here (s re ,f K ) is the unique solution of 

f = exp (/5/i - (3\s/K d ) 

s = Tr[rG K ((3)(l-rG K (f]))- 1 ]. 

Obviously, f K e (0, ||(j k (/?)|| _1 ). Note also that r K and f K satisfy the following conditions 
respectively: 

iTr [r„G„G9)(l - r K GM)~ l ] = ^"^ (2.8) 

^Tr [f„G„G0)(l - fAW)" 1 ] = /3/i ~ 1 A ° grK - (2-9) 

Since by definition (2.5) one obviously gets: G K {(3) ^ G K ((3), the operator G R {(3) 
also belongs to the trace-class Ci(i}). We put the eigenvalues of G K {(3) in the decreasing 
order 

9^ = \\G K {fi)\\>9^>- 
Then, we have gj ^ m (j = 0, 1, 2, • ■ ■ ) by the min-max principle. 

2.2 Approximations of One-Particle Gibbs Semigroups 

Here we establish some relations between Gibbs semigroup {G K (/3)}p>o and the heat 
semigroup {G(/3)}/3>o. Let P K be the orthogonal projection on Sj onto its one-dimensional 
subspace spanned by the vector Q K , and put Q K := I — P K . 

Lemma 2.2 For any r e (0, 1), 

||v / l-e-/[rG K (/3)(l-rG K (/?))- 1 -rG(/3)(l-rG(/3))- 1 ]v /l - e_/ lli - 0,(2.10) 
||a/I - e-fQ K GMQ K {\ - Q K G K (P)Q K )-Wl ~ e~ f ~ #/||i -► (2.11) 

/ioW in £/ie Kmii k — * oo, where \\ ■ ||i stands for the trace norm in Ci(L 2 {R d )). 
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Proof : First, we show the estimates 



\G n K (P;x,y)-G n (P;x,y)\^ 



A' 



I 1 9 I 1 9 

\x\ + \y\ 



Kn d/2-l y K 

TJl -n/3/2n . |„ _ .|2 



if n/?/« ^ 1, (2.12) 
if riP/n ^ 1, (2.13) 



where A' and 5' depend only on d and (3. In fact, by the Mehler's formula one gets for 
ti/3/k ^ 1 the estimate: 



< 



1 



(2vm/3) d / 2 
1 



\G n Mx,y)-G n (P;x,y)\ 

2n(3/K 

\ _ g-2n/3/K 

exp 
1 — exp 



<//2 



exp 



tanh(n/3/2K) ,, l9 i9n 

-(W + M ) 



F-2/1 



2/v 



2k sinh(n/3/fi;) 



(27rn(3) d / 2 
1 



tanh(ri/3/2fi;) .. l9 l9 . 

9 - (g +\y\) 

ZK 

f sinh(n/5 / k) 



exp 



|x-y| 



2k sinh(n/3/fi:) 



< 



< 



(27rn(3) d / 2 

1 
(2vrn/?) d / 2 V k 



V (n^//c) 



x-y 



2k smh(n(3 / k) 



Anf3 \x\ 2 +\y\ 2 n/3 (nf3\ 2 \x — y\ 2 



2k 2k 



k J 2Ksmh(nf3/K] 



exp 



exp 



\x-y\ 



2Ksmh(n/3/K] 

|2 



\x-y\ 



2Ksmh{nj3/ k) 



.4' 



KU 



d/2-1 



\x\ + \y\ 



K 



Here we have used (A. 2) for the first term, (A. 4, A. 6) for the second term and (A. 4, 
A. 8) for the third term at the second inequality and (A. 8) at the last inequality. 
On the other hand for n/3/ k ^ 1, we obtain: 



\G n K ((3;x,y)-n K (x)n K (y)\ 



< 



+ 



1 



{7lK) d / 2 

1 



KK) d l 2 

1 

-exp 



I _ e -2nf3/K 

exp 



d/2 



exp 



tanhf n3 12k) ., l9 l9 . 

-(1*1 +\y\ ) 



\x-y\ 



2 k 



< 



{nK) d / 
1 



tanhf nBI '2 k) ,. l9 l9 . , 
-^T — <\A + \V\ ) ) - exp 

x\ 2 + \y 



I 1 9 I 1 9 

\x\ + \y\ 
2k 



2/tsinh(n/3/fi;) 

\x-y\ 



exp 



2k sinh(n/3//t) 



2 k 



exp 



\x-y\ 



Be- 2n ^ + -(coth(^ 
e V \2k 



2Ksmh(n/3/K) 



\x-y\ 



2fi;sinh(n/3/ 



k)J 



r>i -nj3/2n , |„ _ ,.| 



,d/2 



K 



where we have used (A. 3) for the first term, (A. 5, A. 7) for the second term and (A. 4, A. 9) 
for the third term. 
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Now, let us show the second limit (2.11), notice that the inequality: 

oo 

J2 \G n M x, y) - G n (fc x, y) - Q K (x)QoAy) I 

71=1 

A' ( |x| + M \ 1 



n=l 



2^ ( ^rf/2-1 i 1 ■ K y ' ( w/c )d/2 



+ E (-^72— I 1 

n=\it/P]+l 



B / e -nP/2K . \x-y\ 2 

W 2 I 1 + k ) ^ (2irn(3) d / 2 



<: C'(K 1 - d/2 VK- 1 \o g K)(l + 



I 1 9 I 1 9 

\xf + \y\ 



K 



holds for C, which depends only on d and j3. The integer part is denoted by [■] . Here 
we used the estimates (2.12), (2.13) and 

N (0(1) (d/2>2) 

£^^=1=0(7^/2) (1<)d/2<2 

n=1 [o(logiV) d/2 = 2 

and so on. Now put 

AM : = y/l - e-fQ K GMQ^ ~ Q K GMQ*)~Wl - e~f. 

Then, since \\Q K G K ((3)Q K \\ = e~^l K < 1, one gets in the limit N — > oo the operator- norm 
convergence: 

TV 

A$ = a/1 - e-/ £ Q K G K (/?) n Q« Vl - e-f - A<*> > . (2.14) 

Tl=l 

Recall that Theorem 3.1(i) and Proposition 2.3(i) of [Tib] yield the strong convergence: 

N 

K N = y/l- e-f J2 G(/3) n Vl ~ e-f -> K f > (2.15) 

n=l 

for JV — ► oo. Moreover, we also have the following estimate for the operator norm. 

N 

\\A^-K N \\ = sup |(0, y/l - e-f Y^(Q K G K (P) n Q K - G{p) n )yf\ - e~f 0)| 
< sup / dx I dy\Jl - e~f( x Wl ~ e-fW\(/)(x)<t>(y)\ 

\\<p\\=l J supp f ^supp/ 



N 
X 

71=1 



£ |G™(/3; x, y) - G n (P; x, y) - n K (x)Q K (y) \ 
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^Il\/l-e-/|| 2 sup £|G^C9;a;,i/)-Gr( ) 9; a; ,i/)-n5(a;)no I «(y)H0 

x,2/€supp/ n=1 



for k — > oo uniformly in JV. Here ||yi — e ~'HI stands for the L 2 norm of the function. 
We have used Cauchy-Schwarz inequality at the second inequality. The standard 3e- 
argument yields that A^ — > Kf strongly when k — » oo. 

On the other hand, since the operators A^ K \ Kf are non-negative, we have for k — > oo 
the limit: 

oo 

P^IK-ll^/Hi = TrA^-TrK f = Y,(<Pn,(A iK) -K f )<f )n ) 

n=\ 

oo oo 

= Y. Z}(&»> V 7 ! - e-/(Q K G K (/3)'Q K - G(/J)Vl - e~/0 n ) 

71=1 2=1 

oo oo 

= E Z^"' Vl - e-f{Q K GM l Q* - G{{3) Vl - e-^n) 

(=1 n=l 

oo 

= ^ Tr a/1 - e-/(g K G K (/3) / g K - G(/3) Vl - ^ 

2=1 

oo „ 

= J2 (l-e- f ^)(G l Mx,x)-G l ((3;x,x)-n K (x) 2 )dx^0, 

1=1 "'suPP/ 

where {<j) n }^ =1 is an arbitrary complete ortho-normal system in $j. Note that we can 
exchange the order of summations over n and I, since 

ii.i 

and 

Tr K f = J2(<f>n, y/l-e-fG l (PWl-e-fc(> n ) 

n,l 

are convergent non-negative sequences. Here we used (2.14) and (2.15) in the third 
equality above. Thus, we get hm^^A^ = Kf in €\(Sj) by the Griimm convergence 
theorem, see e.g. [Z]. 

Let us consider the first limit (2.10). Since the identity rG K (/3)(l — rG K (/3)) _1 = 
Y^n=i rn G K (/3) n holds in the operator norm topology and G™(@; x, y) ^ ,(by virtue of 
Lemma 2.2(h) in [Tib]) we get the representation: 

oo 

[rG K (J3)(l - rGM)- 1 ] (*, v) = E ^"^ x ' v)- 

n=l 

Similarly, one gets the representation: 

oo 

[rG{p)[\ - rG{(5))- 1 } (x, y) = ^ r n G n (/3; x, y) . 

n=l 
15 



In fact the series in the right-hand side of the above representations are point-wise 
convergent because the uniform estimates 

1 \d/2 

x(3 txkJ " \2ir/3, 

hold for all k > 0, x, y G M. d and n G N. Here we used (A.l) for the first inequality. 
Hence we obtain the estimate 



/ 1 2 \ d / 2 / 1 \ < 



[rGMii-rGMr 1 ] ix,y)- [rG(/3)(l-rG(/3))-»] {x,y)\ < ^r'RW *, y)-G n (l3; x, y)\ 



n=l 



J^ n -4' (, \X? + \V\\ ^^ 

ra=l n=r«//31+l 



V — 
2np) \tt/3 ttkJ 



A> / |t|2_i_| 9/ |2 x r \K.//3]+l 

^ -A-r d + ' ' l?/l ) + c- 



k(1 — r)\ k / 1 — r 

which tends to zero when k — ► oo, uniformly in x, y G £ for any compact set £, where 
C denotes a constant which depends only on d and on (3. 
Thus, 

sup |[rG re (/3)(l-rG K (/3))- 1 ](x, 2 /)- [rG(/3)(l - rG(/3))- 1 ](x )2 /)| ^0 

x,y£suppf 

holds. We get the first of the announced limits (2.10) by the similar (even simpler) 
argument to the second one. This finishes the proof of the lemma. □ 

2.3 Estimates for the Scaled Mean-Field Interaction 

In the followings, we use the notation B K := 0(k°), which means that there exist two 
numbers c\ ^ c 2 > such that 

Cl n a ^B K ^ c 2 n a . 

Next, we put W K := (G K (/3)) 1/ Vl - e~f and define D K := G K (f3) - G K (f3) = W K W*. 

Lemma 2.3 For large K > 0, the following asymptotics hold: 

(Q K ,D K Q K ) = \\W:Q K \\ 2 = ||v / l-e-/|| 2 (l + (l))/(™) d / 2 , 

llv / r^^n 2 (i + o(i)) 



TrD„ 



d/2 



(tt«(1 - e- 2 ^))' 
and 

g { o K) -g ( o K) = (n K ,(D K -D K Q K \gi K) -QKGMQ^QMn* 



^^(Vl-e-/, [1 + W* K Q K [1 - QkGMQk^QkWJ-Wi - e-f) 

6{K~ d / 2 ) . 
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Proof: For simplicity of notation we suppress everywhere below the index k in g ■ , en , f2( K ) 
and in Q K . 

Note that the /irst equality is a straightforward consequence of definitions (1.8), 
(2.5). The second equality can be derived directly from the Mehler's formula. 

Now by the min-max principle, for d > 2 and k large enough, we obtain from the 
value gi = exp(— /?/«) the following estimates: 

9o = l>g ^(tt,G K ([3)n) = i-(^D K n) (2.16) 

= l-6( K - d ' 2 )>g 1 = l-d(K- 1 )>g 1 . 

Hence the eigenspace of the operator G K (f3) for its largest eigenvalue go is one- dimensional. 
Let Q be the normalized eigenfunction corresponding to go an d let ^ = o^ + 0' , with 
(fi, fi') = 0. Then G K (/3)fi = ^ yields 

aG K (/3)fi + G«(/?)fi' = ag Q + ^ • 

Applying to this relation orthogonal projector P (on Q) and Q = I — P, we obtain: 

ag - a(n, D K Q) - (Q, D K Q') = ag , 
-aQD K n + QG K {f3)QVt' = g tf . 

Since QG K {f3)Q ^ QG K (j3)Q ^ gi < go, the operator g — QG K {(3)Q is positive and 
invertible. It follows from the second identity that 

tf = -a[£ - QGMQl^QDnCl. (2.17) 

Together with the above first identity, 

= (W:^(l-W:Q[go-QG K ({3)Q]- l QW K )W:tt) (2.18) 

follows. 

For brevity, we put 

X' := W* K Q[g - QG K {P)QY l QW K , X := W* K Q[l - QG K {f3)Q}- l QW K 

and 

X := W:Q[g - QG K {f3)QY l QW K . 

Then we get 

X - X' = -XX', 

and hence 

X = X'{l + X'Y l and 1 - X = (1 + X')" 1 . (2.19) 

By definition of W K and (1.8) one gets for large k the asymptotic: 



W*Q = y/l - e-/(7TK)- d / 4 (l + O^ 1 )) . 
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By virtue of (2.18) it implies the representation: 



9o ~ go = (Tr^Vl-e-/, (1 + X')~Wl- e~f)(l + 0(^ 1 )) . (2.20) 

Now, we want to replace in the right hand side of this representation the operator 
X' by X . Note that (2.11) yields \\X - K f \\ ^ \\X - K f \\ x = o(l). Then by 1 - g = 
0(K~ d / 2 ), g -9i= d(« -1 ), l-gi = 0(O and ||W*n|] = 0(«T d / 4 ), we find that 

\\X'-X\\ = {l-~g )\\W* K Q\g -QGMQ}- 1 [l-QGMQ}~ 1 QW K \\ 

<: Ll?±i\ w *qii _ q Gk (p)q]-iq Wk \\ 
go - gi 

which implies \\X'\\ ^ \\K f \\(l + o(l)). Hence (2.20) yields 

(™) d / 2 (so-^o) = (Vl-e-f, (l+XYWl-e-fXl+OiK- 1 )) > lW ^~ *~ ^ (l+o(l)). 

Since the upper bound (7r/t) d / 2 (g — So) ^ (7rft) d / 2 (f2, D re fi) ^ ||a/1 — e~^|| 2 is obvious, 
we get desired estimate go — go = 0(«~ d / 2 ). Notice that the proof of the equality 

so -go = ^wCv^^ 7 , [1 + ^]" 1 v / T^^) 

follows from (2.20) and the estimate: 

\(y/l~e-f, (1 + X'yWl-e-f) - (Vl-e~f, (1 + X)~Wl-e-f)\ 

< || Vi - e-/|| 2 ||(i + X)" 1 !! ||(1 + X')" 1 !! \\X-X'\\= o(l) . 

Here ||vl — e _ ^|| stands for the norm in S) = L 2 (R d ), while the other || ■ || for the 
operator norm on S). This remark finishes the proof of the lemma. □ 

2.4 Evidence of Two Thermodynamic Regimes 

Now we return to the conditions (2.8). We need the behavior of r K and r K to prove the 
main theorem. Here we consider the behavior of r K , which classifies the phase separation. 
That of f K is postponed to Section 4. 

Proposition 2.4 (a) {r K } converges to r* G (0, 1) for k — > oo, where r* is the unique 
solution of 

fj, _ logr« 1 r rjp 

A" #\ /? d i M ^^-r, ' l ' ] 

if and only if f3 d n < ({d)\ . 

(b) K d (l — r K ) — > j3 d \/(P d fi — C(d)X), and hence lim^oorK = 1, if and only if (3 d fi > 

C(d)\. 

(c) lim^^oorK = 1 and n d (l — r K ) — > +oo , if and only if (3 d fi = ((d)\. 
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To this end let us introduce the notation: 



□1 K) : = ^-(n+[0,l) d ) for n G Z d + , 



and define for r G [0,1], v — 1,2 and k G [1,oo), the functions a v { ■ ; r), a!, K ( ■ ; r) on 
[0,oc) d by 



a^0;r) 



re -bli 



(l -re-IPl 1 ) 17 
and by 

4 K) G»;r):= 



if p e a ( K) 

a y (f3n/n; r) if p G D& for ri G Z^ - {0}. 



It is easy to show the following fact. (See [Tib] for detail.) 

Lemma 2.5 There exists a constant c which depends only on d > 2 and (3 > such 
that 

< ai K) (p; r) < 0l/ (cp; l) G ^([0, oo) d ) 

/ioWs /or r G [0, 1] and i/ = 1, 2. 

Remark 2.6 If a series {r K } C [0, 1] converges to r G [0, 1], a!> (-,r K ) — > a^(-,ro) holds 
a.e.. Then the lemma and the dominated convergence theorem yield 

Tr[r K Q K GMQK(l-r K Q K GMQ K )' u } 1 ^ r ^ 






) / ^a { v K) (p;r K )dp-^ (—J / ^a u (p;r )dp 



(3d II V \fl ' kj T \ n^ 

P ' J[0,oo) d V P ' </[0,oo) d 

as k — ► oo. 

Proof of Proposition 2.4: If r K —> r* G [0, 1), then by (2.8) and the above remark one 

gets: 

fi _ logr K 1 /" (k) 



— / aS K) (p;r K )rfp 

P J\0,oo) d 



,d 



lo S r * 1 f , m l f / inj C(d) 

— >—^T + 1m ai(P;n)dp<— a 1 (p-l)dp=— r . 

PA P" ./[o.oo)^ P 1 /[0,oo)« P 

Since equality holds at the limit actually, (2.21) follows. Similarly, if r K — >• 1 and 
K d (l — r K ) — > oo, then one obtains: 

H _ logr K 1 f (k) 



+ ^5rf / aS K) (p;r K )rfp + 

P </[0.oo) d 



A /?A p^p,^ i ^ K> * (l-r K )K 

1 r f u ^ 
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Finally, if n d {l — r K ) — > a > 0, we get 

A* 1 , ((d) . C(<0 



> 



A a (3 d 13 



and a 



(3 d \ 



P d H-((d)\- 



On the other hand, if {r K } does not converge, by compactness we can take two convergent 
subsequences {r Ki } (i = 1,2) with different limits < r^ 1 ' < r^ < 1. Then above 
arguments yield 



l cr r (i) y 



+ lul a 1 (p;r il) )dp= ° g \ + — d / Oi(p;r (2) )dp 



3g 

/5A /3 d J[0,oo) d P^ P~ J[0,oo) d 



which contradicts to the strict monotonicity of the function 

logr 

'[0,oo) d 



h(r) 



(3\ 



— / a x (p;r)dp . 

P J\Q,o6) d 



(2.22) 



Similar arguments are valid for the cases r^ 1 ' = or r^ = 1. If r K —> 1 but K d (l — r K ) 
does not converge, we again get a contradiction. □ 



3 Proof of Theorem 1.2 

3.1 The case \i < fj, c ((3) (normal phase). 

Let us recall that for the weak convergence of random point fields, it is enough to 
prove the convergence of the generating functionals. Therefore, we have to evaluate the 
integral in (2.4). We begin with estimates of the Fredholm determinant in the integrand. 
For all values of x, we have 



Det 


(- 


= 


Det 


= 


Det 



r K G K ((3) 
l-r K G K {p) 

» ^ r K G K (/3) 



J l-r K G K ((3) 

\ , A ■ 2 X r K G K {/3) 
1 + 4 sin — - 



1-fe" 



r K GM 
l-r K G K (J3)J\ 



2(l-r K G K ((3)) 2 
Set a G (d/3,d/2). Then for |x| ^ n~ a , we have 



> 1. 



(3.1) 



Det 
In fact 



1 - r K G K {(3) 



exp 



, xTr r K GM x 2 ^ r K GM 



l-r K G K (/3) 2 (l-r K G K (/3)) 2 J 



(l+0(K d - 3a )) 

(3.2) 



oo 

log of the l.h.s. = - Y^ log (l - (e ix - 1) 



n=0 



(k) 
1 (k) 

1 - r K g\ 
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l)Tr 



r K G K {(3) 
1 - r K G K ((3) 



+ 



-Tr 



Tr 



r*GM 
l-r K G K (p) 

r K G K ((3) 
1 - r K G K ({3)l " \i _ rKg WJ ' 1 - r K G K ((3) " 



where i?i = 0(K d 3a ) since r K — > r* G (0, 1) and 
r*G K {p) 



r <.V% 



+ i?i, 



0(«; d ). 



For the last equality, we recall (2.8) and Proposition 2.4(a). 
Similarly, we have 



l)Tr 



and 



D ix -i \2 



Tr 



r K G K {(3) 
l-r K G K ((3) 

r R G K {{3) 



. _x\ r K G K {{3) 

%X 2) iT l-r K G K (f3) 



0{k 



d—3a\ 



2 -\l-r K G K ((3)J 2 \l-r K G K (P)J 

Thus we get (3.2) and the following lemma. 

Lemma 3.1 

,K d (J3ti,~lo S r K )2/20\fa + Qf K d-3a\\ 



0{k 



d— 3«^ 



-«(/5,^) 



y/1 + /?A K - rf Tr [r K G K {P)(l - r K G K ((3))~ 2 }Det(l - r K G K (f3)) 



Proof : From the above estimates and (2.3), we have 

e -is K x-K d x 2 /2/3\ 

dx ■ 



De t[l - (e- - l)r K G K (/3)0- - r K G K ((3))- 1 } 

^ e - K d z72/3A-is K x+izTr [r K G K (/3)(l-r K G K ^ 



dx exp 



where 



/?, 



/«! , T [A(g u 2 

V/5A + r (l-r K G K (/5)) 2 y 2 

p -is K x-K d cc 2 /2/3A 



;i + o( K d - 3a )) + .R 2 + i?; 



3- 



E|>lt - Q Det[l - (e« - l)r«G K (/3)(l - r K G K (/5))-!] 
for some c > 0, thanks to (3.1), and 



Ofe" 



/?3 



dx exp 



\x\>K~ a 



/5A + r (l-r K G K (/?)) 2 / 2 



(l + 0(K d " 3a )) 



= 0(e- £ ' K ) 

for some c' > 0. Then, the lemma follows from (2.4) and (2.3). 
For S K (/3,/i), we have the following asymptotics: 



□ 
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Lemma 3.2 

S K (/3,/i) 



,K d (/J A t-logf K ) 2 /2/?A/ 1 + 0( w rf-3a^ 



1 + /5A«:- d Tr [f K G K (/3)(l - f K G K (/3))- 2 ]Det(l - f„G K (/?)) 



To show the formula, we note that "tilded" quantities are close to corresponding "un- 
tilded" ones. In fact, the following asymptotics is established. Then it is obvious to get 
Lemma 3.2 by a similar argument. 

Lemma 3.3 For large k one gets: 



(0 

(«■) 

(Hi) 
(iv) 



f K ^ r K , 



r K -r K = 0(K~ d 

-l-r K G K (/?)J 
r K G K (/5) 



Tr 
Tr 



L(l-r K G K (/^ 



Of*; 



Tr 



Tr r f K G K (p) 

h-f R G K {(5) 
r K G K ((3) 



0(k° 



[l-r K G K (f3)) 2 



0(1). 



Proof : (i) Let /i K be the functions on (0, 1) defined by 
h ( ) = l0gr "*" Tr rG<K ^ logr i r 



/3A /t d l-rG K (J3) p\ K d {l-r) 
We also introduce the function /i K on (0, % ) by 

r , x logr 1 rG re (/3) 

M r ) = -^r- — jTr 



-- / a[ K) (p;r)rfp. 



/3\ K d l-rG K (j3)' 

Since G K ((3) ^ G K (j3), h K ^ /i K follows. Obviously /i K and h K are strictly increasing 
continuous functions. Then we have r K ^ r K , because r K and r K are solutions of /i K (r) = 
fi/\ and /i K (r) = /i/A, respectively. (Recall (2.8) and (2.9).) 
(ii) From h K (r K ) = h K (r K ), we have 

f ~ ~ ~ 

— -log— < h K (f K ) - h K (r K ) = h K (r K ) - h R (r K ) 
p\ r K 



K d 



1 



-(GM-GM) 



1 - r K G K (/3) 



0(k 



\-l-r K GM 

Here we have used Lemma 2.3 and the fact that r K is bounded away from 1. The desired 
estimate follows. 

(iii) Since we already know that r K ,r K — > r* 6 (0, 1), we get these estimates readily 
from (2.8) and (2.9). 
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(iv) This can be derived by the telescoping together with (ii) and Lemma 2.3. □ 

Now, let us consider the limit of the ratio S K (/3,/i)/S K (/3, /x) to derive (1.27). From 
Lemma 3.1 and Lemma 3.2, we have 



Ef,, r . [e-Utf] = hm Z Mfi) = lim 



Kr+ooZ K (P,(i) ^Vl + ^lV[r;G^)(l-r;G K (/?))- 2 ] 



Det[l - r K G K (P)] D^t 1 ~ r,G K {(3)} K , 
T>et[l-r K G K {P)]Vet[l-r K G K 



(2/3A)- 1 [(/?M-logr« ; ) 2 -(/3 A1 -logr K ) 2 ] 



Lemma 3.3 yields that the first factor is equal to 1 + 0(k d ). For the second factor, we 
note that 



y/l-e-f 



r«G K (P) 
1 - r K G K ((3) 



\A - e~f - y/l - e~f 



r*G K (P) 
l-r,G K (P) 



y/l - e-f 



l-e-f* 



< 



nGJJ3) n-f K 1 / nG K ((3) 

1-nGM n l-r K G K {p)\l-uG K {p) 



I -e-f 



r* — r 



* ' K 



r*(l -r K ) 



\J\ - e~f 



r,G K {(3) 
1 - r*G K (/3) 



\J\ - e-f 



holds, where || • ||i stands for the trace norm. Together with (2.10), 
Det[l - r K G K (/3)] 



Det[l- f K G K (P)} 



Det[l + r K (G K ((3) - G K ((3))(1 - r K G K {(3)Y 



Det[l + y/l - e-f- 



r K G K ([3) 



1 - f K G K {0) 
follows. Similarly for the third factor, we have 

Det[l - r K G K (f3)} 



1 - e-f] -> Det[l + \J\-e-f 



r*G(P) 



1 - nG{{3) 



'1 - e-f] 



exp 



r K -r K 



Tr 



Det[l -r K G K (f3)} 

r K G K ((3) 
1 - r K G K ((3) 



Det 



l _ r K - r K r K G K ((5) 
r K l-r K G K ([3) 



+0(n- d ) ) = exp (-^^5-(/3/i-logr K )+0(^ d ) 



r K P\ 

where we have used Lemma 3.3 and (2.8). It also follows from Lemma 3.3 that the 
fourth factor is equal to 



exp 



(2/3/i - log r K r K ) log (1 + 



L2/5A 



K ' K 



exp 



r K - r K k" 



K ' K 



r K P\ 



(Ph-\o%t k )+0(k~ 



Thus we get (1.27). Since convergence of generating functionals yields the weak 
convergence of random point fields, Theorem 1.2(i) follows. □ 
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3.2 The case \i > /i c (/3) (condensed phase). 

Let us put: 

^=rvi' ^ =k *( 1+ ^) ail d^=io g (i + ^ 

1 'Ki/j Pj Z Pl 

for j = 0, 1, ■ ■ • . Then it follows from (2.3) and Proposition 2.4(b) that 

f><«>, i +P r^=*^(i + »(.))=*-), 4"»=6 ( o. 

~^ 1 - r K /3 a A 

Since gj = e _/3//re , we have 

pW = 0(«), «$M = d(K^), e[ K) = 0(0 

and 

4 K) < ^ < e[ K) ^ 4 K) < ■ • • . 

Note also that s K = 0(K d ) holds because of (2.3), (2.8). By Remark 2.6 and r K g? e (0, 1) 



Sk = 

j=0 



£p? )2 = o(A 




>g(l >))= ex P(- 
Po 


(k) 

Sk-Po 

(«) 
Po 



(3.3) 



and (A. 10), we get 

oo (k) oo 

n(i-^)=« p(E"«(i-^))=« pf-^+o^ 

j=l Po j=l 

Similarly, we have 



ri( i - p?v (k) - 1)) = e -p ( - e^v w - !) + °( 






exp(-^( SK -^) + 0(^- 2 )) 



,(«)/' „«jW i\ / ^( K )^<5 (,c) 



since pfV - 1) < pfV - 1) = 1/2 for j = 1,2, 
Lemma 3.4 For /arge k we have the asymptotics: 



= ^AZ(, 8 V-CWA)Det(l- r „G K 03)) (1 + ° (1)) - (3 ' 4) 
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Proof : As in the proof of Lemma 3.1, we start with the integral 

-is K x-K d x 2 /2f3X 

X = I dx 



De t[l - (e- - l)r K G K (f3)(l - r K G K (f3))^} 

oo —is K z—K d z 2 /2fi\ / poo poo— i5^ \ poo—iS 1 --^ 

dz rT - = I / - / ) + / = Zi + X 2 . 

Since 1 +Pq ~~ e lz pQ = (1 + p )(l — e 4 (^- 2n7r + 4e o )y the integrand of Xi has simple 
poles at z = 2mr — ie K (nE.li). Let us calculate X\ by means of residues: 

X\ = —2ni 2, R- es 



-is K z-K d z 2 /2(3\ 

; 2mr — ie. 



1 +rf°)(l - e*(-2-+*4 K) )) Yl°° =1 (l - ( e « - l)pf)' 



2vre -^4 K ) + ^4 K ) 2 /2/3A 



— (i + ( c -^/^)) 



27r^Ae- 1+ °^ d > 
~ ^(/3<V - C(d)A) ' 

Here the pole z = —ie^ gives the dominant contribution in the second equality. In the 
third equality, we have used the relations above this Lemma. 
On the other hand, we have 



x 2 |^ 


/•oo e -s K 6 ( -^-K d (x 2 -6 ( - K ^ 2 )/2l3\ 


f-oo^nT=o\^-^ ix+6(K) -M K) )\ 


< \ 


IWX e-o^ 


"\ 


- d ( e *»>- 1)^-1 ° { i} ' 



where we have used 

00 00 

J] |(1 - (e^ W - l)pW)| > ((e* M - 1)# -1) JJ |(1 - ( e * (K) - l)pf )| 

3=0 3=1 

£ ((e 5(K) - 1)^ K) - l)ex P HW( SK -#) + 0(^" 2 )), 
which follows from 

1 - (e"' +5(K) - l)pf = (1+pf )(1 - e^ +5W -4 K) ) 

and <jW -4 K) > 0, <jM -ej° < (j = 1, 2, ■ • • ). Note also that S^p { K) = O^" 1 ) holds. 
From (2.4) and the first equality in (2.3) one gets desired expression for the asymptotics 
oiZ K (/3,u). a 

In order to obtain the corresponding asymptotics for S K (/3,fj), we use the following 
estimates about p* = f K gp/(\ — f K gj) (j — 0, 1, • • ■). 
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Lemma 3.5 

oo oo 

i=i j=i 

Proof : Proposition 2.4(b) , Lemma 3.3(i) and Lemma 2.3 yield 

1 - 6{^ d ) = r K ^ r K < gM- 1 = 1 + 6(«T d / 2 ), (3.5) 

which implies |1 — f K \ = 0(n~ d l 2 ). Note that the argument which shows r K ^ f R in the 
proof of Lemma 3.3 is also valid for the present case. In the variational formula 

9i = sup — - — = sup — — , 

we can use as ip a linear combination of two eigenfunctions K~ d ^(p s (x/ ' \/~k) of G K ((3) 
perpendicular to Q (e.g., with s = (1,0,0, •••) and s = (0, 1, 0, • • •)). Then we get 
gl ^ 1 - O^" 1 ). Together with gf < g[ K) = 1 - d(n~ l ), g[ K) = 1 - 0(k~ 1 ) follows. 

Thus 6{k) = p[ K) >p { 2 K) >■■■ holds. 
Now we get 

oo oo oo 

I E4 K) - E#Va - 1 ] )\ = i 1 - ^ E^V((i - ^f)(! - # } )) 

3=1 3=1 3=1 



= o(^/ 2 )^^V(i-^), 

3=1 

which implies 

oo oo 

j>f =(i+o(« i -^))x;^ ) /(i-yi' s) ). 

3=1 3=1 

On the other hand, because 

oo ~M oo (re) oo (re) ~M _, ,_. 

we have 

., oo ~M -, oo (re) a/ t\ 

K d *-*> 1 _ 5W «; d ^ i _ (") v ; /3 d 

3=1 y 3 3=1 y 3 
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where we recall Remark 2.6. Thus we have n d Y2'jLiPj = P d C(d) + °(1)- Using (2.9), 
we get 

f^ _ lo §^ fi if ~(«) _ ft ((d) 



l (l-fJo K) ) ? X X K *U X P 






which yields the first and the third equality. 

To prove the remaining last bound, it is enough to show that 

P < f ) ^2p < f ) (j = l,2,---) (3.6) 

hold for large enough k, because of Y^jLiPj = 0(n d ). In fact, in the expression 



V 



<«) = r *9j K) l + (f K - r K )/r K 

1 - r K gf> 1 - (f K - r K )g( K) /(l - r K g^ K) 



(f K - r K )/r K = 0(k d/2 ) and \{r K - r K )gf ] /{I - r K gf ] )\ ^ (f K - r K )/(l - r K g[ K) ) = 

0(K 1 ~ d l 2 ) hold. Because of gj ^ g^\ we also have r K (jj /(l — r K gj) ^ pj. Thus we 
get (3.6). □ 

It is obvious now that the next Lemma can be derived along the same line of reasoning 
as the proof of Lemma 3.4. 

Lemma 3.6 For large k one gets the asymptotics: 



2tt{3h pd-l eK *V) li -l OB r K y/20\ 

^(P d fi-ad)X)Det(l-f K G K (P)) 



5 «(/^) = \l /£ (Rdl — r^7^ (l + o(l))- (3-7) 



In order to calculate the limit of H K (/3, /z)/S K (/3, //), we use the following lemma, 
where we put 

gf := (n K , G K (P)Q K ) + r K (Q K , G K {p)Q K (\ - r K Q R G K {[3)Q K )- l Q R G R {p)n R ). 

Lemma 3.7 For large k one gets: 

(i) r K - r K = (1 - #)(1 + o(l)) = O^/ 2 ), 

(ii) l-fJ K) = (l-fJ K) )(l + o(l)). 

Proof: From Lemma 3.5 and Proposition 2.4(b), we have r K % — r K — o(K~ d ). Hence 
(i) follows from f R — r R = f K (l — % ) + t k 9o ~ r Ki ^k — 1 + 0(n d / 2 ) and Lemma 2.3. 
By virtue of Lemma 2.3 we get 

$>-$> = (w* K n K , w:qM k) -QkGMQ.Y 1 - (r- 1 -Q.G K {p)Q K y l }Q K w K w:n K ) 

= (W* K Q K , W:Q K (~g { K) - Q K G K {p)Q K )~\r? - gi K) )(r- x - Q K G K (P)Q R )- l Q R W K W* K n K ), 

(3.8) 

27 



which yields g { K) - g { Q K) ^ since r K < g { K) \ Recall that \\X\\ < 1 and HW^fi^ll = 
0(/s; _d / 4 ). ( See (2.19) and its next line.) Then we also get from (3.8) that 

g { o K) - g ( o K) < j— r^ (w:nw, w:Q K (~ g { K) - Q^MQ^Q^w^) 

1 - r K g{ 

1 _ ~ ~(«) 

^ (^yl^llll^^o II — \. i-r «^o ) U \ K )■ 

1 - »K01 

Hence we obtain the asymptotics (ii): 

1 - ?J K) = 1 - ?J K) + f„(fi#° - g^) = (1 - f K ~g { K) )(l + o(l)). D 

Now, taking into account (3.4) and (3.7), we can find the asymptotics of the gener- 
ating functional (1.25): 

t? r -</,£>! S «(A/i) / K d r K \Det[l-f K Q K G K ((3)Q 

E k r n\e u ' ?/ = „ ' — - = exp ——{2Bu — logr K r K log — 

,AmL J E k (P,h) P V2/3A 1 ^ S k k) & fj Det[l - r K G K ((fi\ 



Det[l - f K Q K G K (/3)Q K ] Det[l - r K Q K G K (/3)Q K ] Det[l - r K G KV 

X Det[l - f K Q K G K (/5)Q K ] Det[l - f K Q K G K (/3)Q K ] Det[l - r K Q K G K ((3)Q n] l ° UJ ' 

(3.9) 
By virtue of Lemma 3.7(i), for the exponent of the first factor, we have 

^-(2/3/1 - logr K r~ K ) log ? = -^(1 - ^ k) )(1 + o(l)) . 
2/3A r K A 

For the second factor, we use the Feshbach formula, which claims 

DetA = Detfi Det(C - K T B- l K), 

where 

B -K\ ( 1 0\ /£ \/l -fi- 1 ^ 



-K T C j == ^-FB- 1 1^ \0 C - K T B- l K) {0 1 

This formula and Lemma 3.7(h) yield 

Det[l -r K Q K GMQ«\ _ 
Det[l - f K G K (J3)] 

1 

1 -f K (^ K) ,G K (/3)^ K) ) - (^\? K G K ((3)Q K {1 -f^GMQ^Q^GM^ 
= l/(l-r^ ) ) = (l + o(l))/(l-f K ^ ) ). 

Detjl -r K GM] 
Det[l -r K Q K G K ((3)Q K ] ^ 
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then Proposition 2.4(b) and Lemma 3.5 yields for the product of factors in (3.9): 

(the 2nd factor) x (the last factor) — ► 1 

in the limit k — > oo. 

Now, since Lemma 3.5 and (2.11) give 

^ ii f, T7 f K Q K G K ({3)Q K r — r — Q K G K (f3)Q K r — 

+IIV1 i-e«G.(/3)g« v /Hl 

* , f . _ n |ivr^7 -j^^- vr^7|i,ii(i - r^GMQ.r\\ 

for ft — ► oo, we obtain the limit: 

Det[l-f K Q K G K (/3)Q K ] 1 



Det[l - f K Q K G K 03)Q K ] Det[l + f K Q K (G K (/3) - G K {/3))Q K {1 - r^GMQ*)- 1 } 

= Det[l + f K yf\ - e-fQ K G K ((3)Q K (l - r K Q K G K (P)Q«)~Wl ~ e^}' 1 -► Det[l + A'/]" 1 

for the third factor in (3.9). Here we have used the cyclicity of the Fredholm determinant. 
Lemma 3.7(i), (3.3), (2.3) and Proposition 2.4(b) yield 

Det[l -r K Q K G K {f3)Q K ] \ 

Det[l - f K Q K G K {f3)Q K ] Det[l - (f K - r K )Q K G K (f3)Q K (l - r K Q K G K (f3)Q „)-*] 

(r K -r K r K Q K G K ((3)Q K ,„ m \ (, _(«) N f d Pn - logr K r K 

= exp Tr +Q(1) = exp (1-^ 



ft- 



r K l-r K Q K G K (/3)Q K ' V V ^V v ™ 'L (3\ l-r„ 

= exp{(l-#)/^(d)(l + (l))//3 d }. 
Thus, by Lemma 2.3 and (2.11) we get for the product in (3.9): 

(the 1st factor) x (the 4th factor) 

^ 2 (l + o(l))p d fi-C(d)\ i 



-0(1' 



exp 



Wl-e-f, (1 + K f )-Wl-e-f)). 



n d/2 pd X 

Now Theorem 1.2(h) follows by collecting in (3.9) the asymptotics of factors that we 
find above. □ 
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4 Proof of Theorem 1.8 

We start with the grand-canonical expectation value of the total number of MF inter- 
acting bosons in the WHT (1.31): 



N K , x (P,n):=K d ff(P,n) 



PE K (P,n)dfj, 



9 S«(/M 



where we use for S K (/?,//) the expression (2.2) after the ^-integration where the values 
of s and r are not fixed yet. The differentiation with respect to \x can be converted 
into differentiation with respect to x in the Fredholm determinant. Then integrating by 
parts we obtain 



N K ,x{P,n) 






where 



R{r, s) 



dx 



xe 



-isx-K d x 2 /2/3\ 



Det[l - (e» - l)rGM(l - rG^/3))- 1 ] 



(4.1) 



dx 



-isx-K d x 2 /2/3X 



De t[l - (e- - l)rG K ((3)(l - rG^/3))- 1 ] 



Then we put s = s K and r = r K such that (2.3) holds, and if we prove R(r K , s K ) = o(l), 
then (1.32) follows as a consequence of this asymptotics. 

To this end notice that for the case (i) /i < /x c (/?), we get R(r K , s K ) = 0(K~ d ^ 2 ) from 
the estimates in the proof of Lemma 3.1. 

For the case (ii) /x > /i c (/3), using the notations of the proof of Lemma 3.4, one 
obtains for the first factor in (4.1): 



dx 



xe 



-is K x-K d x 2 /2/3X 



Det[l 



l)r K G K {(3){l-r K G R {(3))- 1 ] 



-ie X 1 + o(X 1 ) 



Therefore, we get that R{r Kl s K ) = 0(n~ d ). 

The other properties stated in the Theorem 1.8 follow straightforwardly from Section 
2.4 and the line of reasoning developed for the proof of Proposition 2.4. □ 

5 Concluding Remarks and Conjectures 

In the present paper we consider a model of mean-field interacting boson gas in traps 
described by the harmonic potential. For this model we study the position distribution 
of the constituent bosons in the WHT limit by means of the RPF method. 

It is shown that there are two phases distinguished by the boson condensation. In 
one domain of parameters, the resulting generating functional for the RPF is the same 
as for the non-interacting boson gas, for unconventional values of the IBG parameters. 
Whereas in another domain, the generating functional describes divergence of the density 
due to the localization of macroscopic number of particles. 
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Our results are obtained via analysis of the generating functional. We do not intend 
to start the analysis using the characteristic functional here. However, we would like to 
mention a topic on the central limit theorem as a conjecture. 

Let us consider E K ^ :fJi [e~^'^] for small / G Co(R d ) in the sense of the sup-norm || • H^. 
By Theorem 1.2(h) we obtain: 



exp ' 



dx f(x) + hf, (1 + C(/3))(l - GG0))- 1 /) + 

TV I A JM.d ^ 

This leads us to conclusion that the following claim might be true: 

Conjecture 5.1 Let £ K be the random point measure on M, d with distribution given by 
u k,p,h f or K > 0? an d £ be the Lebesgue measure on IR d . Then for k — ► oo the random 
field 



^Kll-fl^iP)/ V^ TT d / 2 X 

converges in distribution to the Gaussian random field on M. d with covariance 

(l + G(/3))(l-G(/3))- 1 . 

We finish by some remarks about the method of the RPF approach to the BEC in 
the WHT limit used in the present paper: 

(i) It could be applied to a general "non-quadratic" mean-field interaction [/$ := 
k <&(n/K ), where $:xGMi— > 1 is a piece-wise differentiable continuous function 
bounded from below [TZ] , as well as to the van der Waals particle interaction, which is 
more local than the mean-field [deS-Z]. We guess that for this kind of interaction the 
particle distribution will spread as n a with some large a even for condensed particles, 
(ii) The method has to be compared with the scaled external field perturbation of BEC 
considered in [deS-Z] , [Pu] . We suppose that this could clarify the concept of the choice 
of "effective" volume, since it is important for description of the local particle density 
measured in the WHT limit BEC experiments as well as its definition of the mean-field 
interaction [DGPS], [LSSY], [PeSm]. 
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A Miscellaneous formulae 



X 



1-e- 



X 



^l+x^2V(2x) for x > 0. (A.l) 



2x N d/2 



1 - e" 2 * 



< ( :z ^ ) - 1 ^ Ax for x E [0, 1], (A.2) 



1 



d/2 



^ - r - 1 < Ee - ^ for x ^ 1, (A.3) 



1-e- 5 
where A and .B are constants depending only on d 



1 - e~ x ^ x for x ^ 0. (A.4) 

e" x - e~ y ^ - ( - - 1 J for < x ^ y. (A.5) 



tanhx ^ x A 1 for x ^ 0. (A. 6) 

2e 
e-1 



2e 
cothx - 1| < e _a; for x > 1. (A.7) 



sinhx 9 . , , . 

1 ^ ^ 1 + x 2 for x G [0, 1]. (A.8) 

-\~ < : 2 , , 2 e" x for x^l. (A.9) 
smnx 1 — 1/e^ 



x 2 



log(l - x) - x| ^ — for x < 1. (A. 10) 

Zl 1 Ju J 
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